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Abstract 

We present here a simple model of radiative gravitational collapse with 
radial heat flux which describes qualitatively the stages close to the formation 
of a superdense cold star. Starting with a static general solution for a cold 
star, the model can generate solutions for the earlier evolutionary stages. The 
temporal evolution of the model is specified by solving the junction conditions 
appropriate for radiating gravitational collapse. 
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1 INTRODUCTION 



The theory of gravitional collapse has many interesting applications in astrophysics 
where the formation of cold compact stellar objects are usually preceded by a period 
of radiative collapse. The problem of gravitional collapse was first investigated by Op- 
penheimer and Snyder [|IJ , who considered the contraction of a spherically symmetric 
dust cloud. Here the exterior spacetime was described by the Schwarzschild metric 
and the interior space time was represented by a Friedman-like solution. Thereafter, 
Vaidya derived the line element which describes the exterior gravitational field of 
a spherically symmetric radiating mass. It then became possible to model the inte- 
rior of radiating stars by matching the interior solution to the exterior space time 
given by Vaidya []3], f|, [|, ^ [7J]. The junction conditions for a spherically symmetric 
shear-free radiating star was completely derived by Santos ||. The crucial result 
that follows from Santos is that the pressure on the boundary of a radiating sphere 
is nonvanishing in general. Subsequently, many models of radiative gravitational 
collapse with heat flow were found by utilising these junction conditions. (See, e.g. 
0.) In particular, special attention was given to models in which an initial static 
stellar configuration began collapsing by dissipating energy in the form of a radial 



heat flux fL0| . The initial static configuration was taken to be an exact solution of 
the Einstein field equations. 

The aim of this paper is to consider the evolution of a star undergoing radiative 
gravitional collapse with its final state being that of a superdense star. We wish 
to study this evolutionary process starting from the final nonradiating state and 
interpolating to earlier times when it was emitting radiation. 

In Section 2 we present the relevant background material and the line element 
for the interior space time. The Einstein equations for an energy momentum tensor 
with heat flux are solved without assuming a particular form for the final static 
configuration. The Vaidya solution is introduced in Section 3 and the matching 
condition is utilised to determine approximately the temporal evolution of the model. 
In Section 4 we use the solution of Mukherjee et al (Tl| to describe the final state 
of the star and find out its evolution by making use of the final static configuration. 
Although we are using a simple model, the results obtained analytically are expected 
to provide the general trends of results for realistic stars. In Section 5, we summarize 
our results. 
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2 INTERIOR SPACETIME 



Let us assume that the interior spacetime of a star is represented by a line element 
of the form 



ds 2 



-A 2 {t,r)dt 2 + B 2 {t,r){dr 2 



(1) 



The choice of a shear-free metric is motivated by the simplicity of the resulting 
equations, which remain tractable throughout. Moreover, as observed by Bonnor 
et al jnj, it is possible to show through Raychaudhuri's equation that the slowest 
possible collapse is for shear-free fluids. Thus our results will be appropriate in any 
case if the collapse is not very fast. Also, the choice of a shear-free metric will later 
on help us to apply the same junction conditions derived by Santos || in our model. 
The energy momentum tensor for the interior matter distribution is taken to be 



T^ v = (p + p)u^u u + pg^ + q^u u + q u u h 



(2) 



where the flow vector satisfies 



0. 



vJ 1 being a timelike four velocity vector. 

Following Bonnor et al ||10||, we choose the metric functions as follows: 



A(r,t) = A {r) + ea(r)T(t) 
B(r,t) = B (r) + eb(r)T(t) 

and the energy density p and the isotropic pressure p as 

p(r,t) = p (r) + ep(r,t) 
p(r,t) =p (r) +ep(r,t) 



(3) 
(4) 



(5) 
(6) 



The radial heat flux is of the order of e (0 < e 1). However, unlike the case 
studied by Bonner et al Aq and Bq describe here the final static solutions of the 
cold star. 

Einstein's field equations for the static configuration give the relations: 



i 



2 { ^ 



(7) 
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(8) 



where prime (/) denotes differentiation with respect to r. 
The pressure isotropy equation is given by 



' \ 2 



(9) 



We assume that the presure isotropy condition is satisfied for a known static so- 
lution (A ,B ). The perturbed field equations up to first order in e can be written 



as 



P 



P = "3poiT + ^ [_ (| ) 2 b + 2 (| _ |) V _ 2 ^ 



(10) 



T " 2^(11) 



9 — H? 



(12) 



where an overhead dot denotes differentiation with respect to t. 

The condition of pressure isotropy for the perturbed matter distribution can be 
written as 



(*)' + (*)T-'[(*)'+(*)1(£ + t) 

(*)' + (*)' 



_|_ AllQ. _a_ 



(13) 



The general solution to this equation is 

/ f-2_V A 2 \ ' 



Bo 



or, equivalently, 



drdr 



(14) 



(15) 



where Co and Ci are constants. 

Equation (14) or (15) gives a relation between a(r) and b(r). For specific cal- 
culations, however, it will be helpful to follow an alternative approach. We write 

= X(r) and (J^) = Y(r) and rewrite equation (13) as 



= 0. 



c a x+y 

Equation (16) can be integrated easily if we assume 

X a dg(A ) 

X+Y ^° dA 



(16) 



(17) 



which gives 



X = hrAlB^ye- 4 ^ (18) 

and 

Y ={^-^) X (19) 

where k\ is an integration constant and y = dg j^ ■ 
Thus we obtain 

% = hj rAlBlye~ i9 ^dr + k 2 (20) 

where k 2 is another integration constant. 

The right hand side of equation (20) can be integrated for different choices of the 
function g(A ). It is instructive to consider a simple case: 

g(A ) = ±\nA (21) 

This gives 

* = i = ^/^r + fe. (22) 

It may be useful to calculate the total energy entrapped within the surface S of 
the star. Up to first order in e, this is given by 

m(r-E, t) = m (r s ) + em(r s , t) (23) 

where 

mo ( rs) = _( r 2 jB / +r 3^ (24) 



mi 



M^([-r%'-r%(2^-l)]T(t)) E (25) 

where rs corresponds to the boundary. The evaluation of mo(rs) and m(rs,t) will 
be pursued in Section 4. 
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3 JUNCTION CONDITIONS 



The boundary of the collapsing star divides spacetime into two distinct regions, the 
interior spacetime described by the metric (1) and the exterior spacetime. Since the 
collapsing star is radiating energy, the exterior spacetime is described by Vaidya's 
outgoing metric 



ds 2 



(» 



2m(v) 



dv 2 - 2dvdr + r 2 (d9 2 + sin 2 



(26) 



m(v) being an arbitrary function of the retarded time v. The solution (26) is the 
unique spherically symmetric solution of the Einstein field equations for radiation in 
the form of a null fluid. The Vaidya solution is often used to describe the exterior 
gravitational field of a radiating star |3], |], |5], |6|, [12[ . 



In order to facilitate the smooth matching of the interior space-time to the Vaidya 
exterior, we utilise the junction conditions deduced by Santos |§. We can rewrite 
equations (11) and (12) as 



P 



-2PoiT + 2^(aT-T 



where 



a 



Abe 

A'2 R2 



4 2 

bB 



(3T 



B, 



By 



A 



+ ( _L I 



A'2 

26 



{a b ) ■ 



(27) 
(28) 



(29) 



(30) 



To find T(t) we make use of the junction condition p s = (qB)^ together with 
(Po)s = which gives 

a E T - f = 2/feT (31) 

where E represents the boundary of the star. The solution of (31) is given by 



T(t) = T exp 



which satisfies the boundary conditions 
T(t)| t=oo = 0and T(t)\ t=0 



(32) 



(33) 



where T is a constant. Since we expect T(t) to decrease as t increases, we have 
«s > 0. (see later.) 

6 



4 FINAL STATIC SOLUTION 



We assume that the line element (1) represents a static solution at large times i.e. 
as t — > oo, and the metric coefficients are then represented by A an d B , the static 
part of equations (3) and (4), respectively. 



For the static solution we take the solutions of Mukherjee et al [1 1] of the Vaidya- 

(34) 



Tikekar [|T3| model for a class of static spherically symmetric stars, viz. 



- e Mr) dt 2 + e 2M(f)d f 2 + f 2(^2 + ^2 







where 



and 



oMf) 



1+Xf 2 /R 2 
l-f 2 /R 2 



Mr) = D 



cos[(n+l)C+H 
n+1 



cos[(n— 
n— 1 



(35) 



(36) 



In equations (35) and (36), D, u and R are constants, C = cos 1 z, n 2 = A + 2 and 

^ 2 = (^r)(i-ft)- 

The solution satisfies strong and weak energy conditions as well as the causality 



condition and is valid for any value of A > A [|Tl|]. The boundary of the star b is 



determined by imposing the condition p = at r = b. At the boundary of the star, 
the solution is matched to the Schwarzchild exterior solution. The model has four 
parameters A, R, D and oj, of which three, say R, D and uj, are determined by the 

matching condition (p = at f = b and the continuity of two metric coefficients). 
The choice of the remaining parameter A, then, fixes the relevant equation of state 
of the star. Thus, in the present model, for a given mass and radius, there is a one- 



parameter class of equations of state parametrized by A. Previous works [0, [15|] show 
that this model can describe cold compact stars like Her X-l and SAX J1808. 4-3658. 
The model also has a scaling property which allows a description of a class of stars 



with the same compactness [16|. Moreover, the stability of the stellar configurations 



for various A can be verified easily in this model, as shown in fL5| . Thus although 
the solution provides a toy model of a star, its simple analytic form is found to be 
very useful in studying the gross features of a star with an equation of state specified 
by the value of A. 

However, to make use of this solution, we need to transform the metric (1) to the 
form given by (34), i.e. we need to change the radial coordinate r to f. 
A comparison of equations (1) and (34) yields 



r = rBo(r) and e^dr = Bo(r)dr (37) 
Using equations (35) and (37), we obtain the inverse transformation relation as 



r = & 3 exp 



^■ 1 fo&)-^* rl G/S£ - »*) 



(3f 



where x — cos 1 (^) an d ^3 is an integration constant. 

The integration in equation (22) is over r, but using equation (37) we can integrate 
it over f so that the final results are expressed in terms of f. We obtain 

JL = Jl = _J|_#2^T Jsin-^/x) + ZaVl - ^ 2 ] + h (39) 



where/ = V / A+T' X = V 1 "^ • 

To get a feel for the behaviour of the solution, we need to consider specific numeric 
cases. As an example, we consider a star whose final state has a mass m = O.88M 
and radius b = 7.7 km. Note that, these values fall well within the estimated mass 



and radius of the well known X-ray pulsar Her X-l ]T7|]. Moreover, with a particular 
choice of the parameter A (A = 100), it has been shown elsewhere JH] that the 
equation of state obtained using the model discussed in Section 4, agrees accurately 
with the equation of state obtained by Horvath and Pacheco [18[] for a quark-diquark 



mixture. Matching the static solution to the Schwarzchild exterior solution and 
using the boundary condition po = at f = b (note that in the evolutionary stage 
pressure does not vanish at the boundary), we calculate the values of the constants 
as R = 108.779km, u = 2.58577 and D = 26.6039. 

We, then, choose the constants ki, k 2 & k 3 in such a way that our model describes 
the expected early evolutionary stage of the star. The time dependence of T(t), heat 
flux q(t) and mass m(r,t) are shown in figures 1, 2 and 3, respectively, where we 
considered two exemplary cases: (1) k\ = 1, k 2 = 10 5 and k- A — 1 (solid lines) and 
(2) k\ — 1, k 2 — 10 5 ' 2 and k 3 = 1 (dashed lines). Arbitrary choices of these values 
may not give realistic results. In figures 1 and 2, we find that both T(t) and q(t) 
decreases with time, as expected. In figure 3, the mass of the star decreases with 
increasing time and it can be shown that as time goes to infinity the mass saturates 
to its final static value of 0.88 M . 



8 



5 CONCLUSIONS 



Our model gives a description of the evolution of a radiating star. For a known static 
configuration of a star, this model generates solutions for the earlier stages of the 
star by a perturbative approach. This may be looked upon as complementary to the 
approach of Bonner et al [|l(| . Possibly, the two methods, when combined carefully, 
will be able to give a total picture of the radiative collapse of a star, which ends up as 
a cold compact star. In our method, we have made use of the general solution given 
by Mukherjee et al [O]. It has been shown elsewhere that the static solution can 



describe very compact cold stars, with an equation of state relevant for deconfmed 
quarks [15] as well as a mixture of quarks and diquarks |14]. The possibility that a 
compact star could have a deconfmed quark phase in the core, surrounded by a crust 
of baryons and separated by a layer where the relevant density permits a hadronic 



phase transition, has also been considered |L9| in the static model. Since a variety of 
final states are now available in the static model, the approach presented here may 
be useful. 
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Figure 1: Plot of 4^ against t. (h = 1, k 2 = 10 5 & k 3 = 1 (solid lines); h = 1, 
k 2 = 10 5 - 2 L k 3 = l (dashed lines).) 
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Figure 2: Variation of heat flux q as a function of t. (ki — 1, k 2 — 10 5 & k 3 = 1 
(solid lines); fci = 1, k 2 = 10 5 2 & k 3 = 1 (dashed lines).) 
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Figure 3: Evolution of mass of a star whose final static configuration has mass 
m = 0.88 M . (£4 = 1, k 2 = 10 5 & k 3 = 1 (solid lines); ^ = 1, A; 2 = 10 5 - 2 & k 3 = 1 
(dashed lines). ) 
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